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ABSTRACT 

We present here a detailed, self-contained treatment of the mathematical formalism 
for describing the theory of polarized anisotropy in the cosmic microwave background. 
This didactic review is aimed at researchers who are new to the field. We first develop 
the mathematical tools for describing polarized scattering of CMB photons. Then we 
take the reader through a detailed derivation of the line-of-sight formalism, explaining 
the calculation of both temperature and polarization power spectra due to the scalar 
and tensor perturbations in a flat Universe. 

Subject headings: cosmology: theory - cosmic microwave background - polarization 



1. Introduction and Brief History of CMB Polarization 

In the gravitational instability paradigm of structure formation, density perturbations in the 
early Universe grow into the large scale structure we observe today. The presence of the pertur- 
bations at the epoch of recombination causes the Cosmic Microwave Background (CMB) to be 
polarized through Thomson scattering. The degree of polarization in the CMB was predicted to 
be weak (about 10% or less than the temperature anisotropy at small angular scale, even smaller 
at large angular scales, see e.g., Kosowsky 1996), the recent detection of polarization in the CMB 
by the Degree Angular Scale Interferometer (DASI, Kovac et al. 2002, see also Leitch et al. 2005) 
and the measurement of the large-angle power spectrum of correlations between temperature and 
polarization anisotropics by Wilkinson Microwave Anisotropy Probe (WMAP, Bennett et al. 2003; 
Kogut et al. 2003) were important confirmations of the paradigm (e.g. Hu & Dodelson 2002, and 
references therein). 

There are several other reasons why studies of CMB polarization are important. First, al- 
though studies of the temperature anisotropics in the CMB have provided strong constraints on 
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many of the cosmological parameters, degeneracies between some combinations of parameters exist 
(e.g., Seljak 1997). In principle, observation of polarization can help break some of the degenera- 
cies, such as that between the overall amplitude of the temperature power spectrum and the epoch 
and degree of reionization (Hu & White 1997a). Second, some physical mechanisms (e.g., gravity 
waves) only contribute to very large angular scales, where our ability of extracting cosmological in- 
formation is limited by cosmic variance. The additional information in the polarization anisotropies 
can add valuable information for studies of cosmological physics on super-large scales (e.g., Zal- 
darriaga & Seljak 1997). Related to this point, we note that different perturbations modes (scalar, 
vector and tensor) give distinguishable polarization patterns (Hu &: White 1997a). Thus, from 
the construction of polarization power spectra we ought to be able to investigate the nature and 
origin of perturbations presented in the early Universe. Finally, unlike temperature anisotropies 
which are affected by various physical effects that occur between the last scattering surface and 
present, polarization provides us with direct probe of the last scattering surface. It therefore helps 
distinguish the contributions to the temperature power spectrum from gravitational potential and 
peculiar velocities (e.g., Kamionkowski, Kosowsky, &; Stebbins 1997). 

Although the importance of the CMB polarization was recognized very early on (e.g. Rees 
1968; Caderni ct al. 1978; Negropontc & Silk 1980; Kaiser 1983), formalisms that enable fast and 
accurate computations of the polarization field on the whole sky were not realized until late 1990s. 
One main reason is that, in applying the Stokes Q and U parameters to describe the polarization 
field, a fixed coordinate frame is needed. If one models polarization maps on the whole sphere, this 
requirement complicates the computation. One way to avoid this is to expand the polarization field 
by tensor harmonics or their relatives, the spin-weighted harmonics. This was first done by two 
groups (Zaldarriaga & Seljak 1997; Kamionkowski, Kosowsky, & Stebbins 1997). In this review we 
follow the formalism developed by Zaldarriaga &: Seljak (1997, hereafter ZS97) to reduce the "spin 
weight" of the polarization fields to zero using the so-called spin raising and lowering operators. 
In this framework the calculation of the polarization power spectra becomes greatly simplified and 
very similar to that of the temperature power spectrum. The method applies to all angular scales, 
and is therefore applicable to the analysis of all-sky surveys. Due to its calculational convenience 
the spin weight formalism is the technique of choice for codes solving the linearized cosmological 
Boltzmann equations such as CMBFAST or CAME ^ 

In this review we first provide a self-contained overview on the basic physics of polarization 
and their mathematical description. Then we detail the derivations of several important results 
of ZS97 (presented in their sections II to IV). We mainly provide the mathematical details of 
the calculations that lead to temperature and polarization power spectra. This is intended to 
help build the mathematical toolkit for new workers in CMB theory. We assume some basic 
familiarity with electromagnetism. One subsection introduces polarization in a quantum mechanical 
notation, however, this is not required in order to understand what follows. Tutorials on the physical 
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interpretation and visualization of CMB polarization have been available for some time and are 
summarized in § 5. 

This review is organized as follows: in § 2 we present some basic ingredients needed to under- 
stand the generation of polarization at the last scattering surface, including the Stokes parameters, 
the spin-weighted spherical harmonics, and the Thomson scattering. In § 3 wc introduce some basic 
definitions for calculations of temperature and polarization power spectra, then formally calculate 
the power spectra induced by scalar and tensor perturbations. We briefly summarize our results in 
§ 4 and provide references for further reading in § 5. In Appendix A we discuss the parity of the 
polarization modes E and B; in Appendix B we describe the line-of-sight integral solutions to the 
Boltzmann equation. 



In this section we discuss three topics which are helpful in understanding the physics of CMB 
polarization: the Stokes parameters, the spin-weighted spherical harmonics, and Thomson scatter- 
ing. In the discussion of the Stokes parameters, we present both a classical (§ 2.1.1) and quantum 
mechanical (§ 2.1.2) picture for describing electromagnetic radiation, then consider observations 
with finite bandwidth of radiation from uncoherent sources (§ 2.1.3). In our section of spin- 
weighted spherical harmonics, we simply state the basic features of this family of functions without 
proof. Finally, the discussion on Thomson scattering focuses first on the scattering matrix, then 
on the generation of polarization at the last scattering surface. 

We follow closely the discussions in Shu (1991), Chandrasekhar (1960), and Rybicki & Light- 
man (1979) in sections 2.1.1, 2.1.3 and 2.3.1; the material presented in sections 2.1.2 and 2.3.2 
follows the discussions in Kosowsky (1996, 1999); § 2.2 is based on the appendix of ZS97, which 
in turn is drawn from Newman &; Penrose (1966) and Goldberg et al. (1966). Note the difference 
in the sign convention for rotation between that of ZS97 and the original reference; it is chosen so 
to match the practical convention used in the astronomical literature. As far as possible we try to 
stick to a single convention in this review, pointing out where our references differ along the way. 



Consider a plane electromagnetic wave propagating along the z direction. Its Fourier decom- 
position can be expressed as 



2. 



Polarization Basics 



2.1. 



Stokes Parameters 



2.1.1. Classical Description 
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where £ and ^ are real quantities denoting the amplitudes and phases in the two transverse direc- 
tions marked by the unit vectors x and y. The angular frequency of the wave is a; = kc. The real 
part of a given k mode is 

Efe = x£x cos{kz — u!t + (f)x) + y£y cos{kz — u!t + (f)y). (1) 

On the X — y plane, the tip of the electric vector will trace out an ellipse as a function of time. 
Let the angle between the major axis of the ellipse {x') and x-axis be x, i.e., 




cos X sm X 
- sin X cos X 




(2) 



We can choose the zero point of time so that Ejt(t = 0) is purely in the x' direction: 

Ek(t = 0) = x' El cos ujt + y'E2 sin ut, (3) 

where 

E! + Ei=£'^+£'^^£i (4) 

and 

Ei = £o cos /3, E2 = £0 sin/3, (-7r/2 < /3 < 7r/2). (5) 

The ellipticity angle j3 determines the shape of the ellipse. For example, for j3 = ±7r/4, the ellipse 
becomes a circle, and the wave is circularly polarized; if /3 = 0, ±7r/2, the ellipse "collapses" into a 
line, and the wave is linearly polarized (also see below). 

The quantities in Eqns (1) and (3) can be related with the help of Eqns (2) and (5); the 
coefficients of xcoswf give £xCOS(t>x = E'lcosx = £^0 cos /3 cos x, those of xsinwi give f^jSin^j. = 
— £^2sinx = — £"0 sin/3sinx, those of y cos tot give £y cos (f)y = Eisinx = £^0 cos /3 sin x, those of 
y sin cot give £y sin (f)y = E2 cos x = ^0 sin (3 cos x- 

The Stokes parameters are then defined as follows: 

/ = £l + £l = £l (6) 

Q = £l-£l = £lcos2f3cos2x, (7) 

U = 2£x£y cos{(f>y — <j)x) = £0 cos 2/3 sin2x, (8) 

V = 2£x£ysm{(t>y-4>x)=£lsui2(3. (9) 

Note that when the wave is monochromatic these parameters arc related by = Q'^ -\- U"^ + V"^ . 

The Stokes parameters defined in this way are all real quantities. The / parameter measures 
the radiation intensity and is positive. The other parameters describe the polarization state and 
can take either positive or negative sign. The parameters Q and U measure the orientation of the 
ellipse relative to the x-axis and define the polarization angle 

X-^tan-^ (10) 
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and the polarization vector 

P = (Q2 + [/2)l/2^/. (11) 

A few comments on the various sign conventions of Stokes Q and U are in order. In the usual 
spherical coordinates, Q > for a N-S (longitudinal, ieg) polarization, Q < for a E-W (azimuthal, 
±e<j(,) component; on the other hand, a NE-SW (±(e0 — e6i)/\/2) polarization component means a 
positive U while a NW-SE (±(e0 + ee»)/\/2) component means a negative one (Hu &: White 1997a). 
The last parameter, V, measures the relative strengths of two polarization states: linear polarized 
light has F = 0; light for which Q = U = Q and F > (F < 0) is right- (left-) circularly polarized. 
In particular, for unpolarized, or "natural" light, Q = U = V = Q. 

It is apparent that, with respect to a rotation about the z axis, /, V and = + U"^ are 
invariant, because they are independent of %, the angle between the polarization vectors {x' Sz y') 
and the artificially chosen coordinates {x &l y). However, Q and U transform in the following way 
under a counterclockwise rotation of the x — y plane through an angle a about the z— axis, as can 
be obtained by letting x ^ x — a in Eqns (7) and (8): 

Q' = Q cos 2a + U sin 2q;, 

U' = -Qsm2a + Ucos2a. (12) 

A compact way of writing these is through a combination of Q and U: {Q ± iU)' = exp(^2iQ!)(Q zt 
iU). In fact this is the property that makes {Q±iU) a spin-ib2 quantity, a fact that we will exploit 
heavily in section 2.2. 



2.1.2. Quantum Mechanical Description 

A quantum mechanical description of the Stokes parameters^, which is closely related to the 
photon density matrix, can be found in Kosowsky (1996). Consider two orthonormal linear polar- 
ization basis kets |ei) and |e2). An arbitrary state vector can be spanned by these basis kets 

\e) = aie'^'\ei)+a2e'^^\e2). (13) 

The Stokes parameters can be defined as the expectation values of four operators 

-f = |ei)(ei| + |e2)(e2|, 

Q = |ei)(ei| - |e2)(e2|, 

U = |ei)(e2| + |e2)(ei|, 

V = z|e2)(ei| -z|ei)(e2|. (14) 



detailed understanding of this subsection is not required to understand the remainder of this review - the 
reader unfamiliar with Dirac notation can skip to 2.1.3 for the first reading. 
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For example, the expectation value of / is 

I = {i) = (e|i|€) = [ale-'^'{ei\ + a^e-^^^(e2|) / {aie'^'\ei) + a2e'^^\€2)) = \ai\^ + |a2|^ 

consistent with the definition in Eqn (6). 

Since polarization represents a mixture of the two degrees of freedom of a photon, the photon 
density matrix has the components(e.g. Sakurai 1985) 

p=f^" ^-y (15) 

y P21 P22 J 

The components of p can be expressed in terms of expectation vahics of the operators defined in 
Eqn (14). We use the definition of operator expectation values in terms of the density matrix; 
e.g. Q = {Q) = ii{pQ) = pii - P22- Similarly, 7 = + ^22, U = {U) = pi2 + P21, and 
V = {V) = i{pi2 — P2i)- Prom these we can express the density matrix in terms of the Stokes 
parameters and the Pauli matrices af 

1/ I + Q U-iv\ 
2yu + iV I-Q ) 

^ {II + Uai + Va2 + Qc73) , (16) 

The density matrix defined here is equivalent to the intensity tensor defined in ZS97: the differ- 
ence is that the intensity tensor is described in terms of fractional change in the CMB temperature 
(recall that the intensity is proportional to the fourth power of temperature, I oc T^). 



P = 



where 1 is the 2x2 unit matrix. 



2.1.3. Practical Description 

So far we have only considered monochromatic waves. In reality, measurements of electromag- 
netic waves usually actually measure averages over several oscillation periods and over a range of 
frequencies. We thus define the averaged Stokes parameters as 

Q = {£l-£l)= I cos2(3cos2x. 

U = 2 {£x£y) cos{(l)y — ^x) = I COS 2/? sin 2%, 

V = 2{£x£y)sm{(t)y-(l)x) = Ism2^. (17) 



Note that the time scale for variations in £i{t) and (pi{t) {i = x,y) must be large compared to the 
period of the wave. 
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An important property of the Stoles parameter is that they are additive for incoherent su- 
perpositions of waves. Light from astronomical sources is therefore not expected to be completely 
elliptically polarized, for they generally come from different regions of the source, with different 
polarization amplitude, vector, and phases. Suppose we receive a beam composed of a mixture 
of independent, elliptically polarized light, = J2j^k^ where different components do not 
possess coherent phases with each other. Upon averaging over time and bandwidth, we have 

j j j j 

2 2 2 2 

In this case, it can be shown that / > Q + U + V (Chandrasekhar 1960, p. 32; Rybicki & 
Lightman 1979, p.67). Because of this property, it is always possible to decompose the observed 
radiation (7, Q, U, V) into two components: one completely unpolarized, the other elliptically po- 
larized. Specifically, for the polarized part, (Jp, Qp, Up, Vp) = {{Q^ + C/^ -|- y^)^/^, Q, [/, V); for the 
unpolarized part, (/„, J7„, Vy) = {I — Ip-, 0, 0, 0). Finally, the fractional polarization is given by 
the ratio Ip/I- In what follows we will only focus on the polarized part of the radiation, and ignore 
the subscript p. 



2.2. Spin— Weighted Spherical Harmonics 

Having described the properties of polarized light, we now develop the mathematical machinery 
necessary to represent angular distribution of the polarization of the CMB on the celestial sphere. 
The representation that is most convenient for making contact with cosmological theory is in terms 
of spin- weighted spherical harmonics. 

A quantity r] that transforms as r/' = eyi'p{— s ia)ri under a rotation of an angle a is defined as 
of spin-weight s. Examples we have already encountered are: the scalar temperature anisotropy, 
which is a spin-0 quantity, and Q±iU which is a spin-zb2 quantity. It is important to note that the 
rotation here does not mean a global rotation that changes labeling of coordinates of the sphere 
(the sky). What is meant is simply a rotation in the plane tangent to the point of interest on the 
sphere, ie., a rotation of the coordinates defined on that tangent plane. 

Just as it is the case that a scalar (spin-0) function on a sphere can be expanded into a series of 
spherical harmonics Yi^{0,<j)) (also spin-0 functions), a spin-s function can be expanded in spin-s 
spherical harmonics sYim{Q ■, <i>) ■ For each |s| < Z, they form a complete, orthogonal basis on the 
sphere: 

j dn,Yir^,{e, (P)sYUe, 0) = Sm s^'m, (is) 

sYCm{e, ct>)sYlmi0', </-') = S{cf> - cf>') d{cOS 9 - COS O'). (19) 

l,m 

Spin-weighted spherical harmonics were devised by Newman k, Penrose (1966). They have 
been applied mostly in theories of multipole expansion of gravitational waves. In this section we 
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foUow the discussion in ZS97, simply state some of their properties without proof. See Newman & 
Penrose (1966); Goldberg et al. (1966); Penrose & Rindler (1986); Hu & White (1997b) for detailed 
descriptions of the function. Also sec Thorne (1980); Kamionkowski, Kosowsky, & Stebbins (1997) 
for relations between this family of functions with other tensor spherical harmonics. 

There exist a pair of operators, 9 and 0, known as the spin raising and lowering operators, 
respectively (in this review, for notational convenience, they are denoted as: () = S, t> = 6). These 
operators have the property of raising/lowering the spin weight of a function; denoting a quantity 
in a frame rotated ip from the original frame by a prime, we have (Us/)' = exp(— i(s + 1)'0) (tts/)) 
and (bs/)' = exp(— z(s — l)'tp) (b^/). Their explicit form is 



UfiO,4>) = -sin^ 

Kf{e,ct)) = -sin- 



sm ( 



df, 



sm( 



-A, 



sm-'esfie,4>), 

sin' esf{0,ct>)- 



Suppose we have spin-ib2 functions ±2/(^5 0) whose (/)-dependence satisfies d^gf 
Acting twice with the spin raising/lowering operators on ±2/ gives (ZS97) 

2 



b'2/(/^,</') 



1 — 



[(l-/x2)2/(/X,0)], 



f_2/(/x,</.) = (-d^-^^y [{l-^^^)_2f{^i,c|>)], 



(20) 
(21) 

irusf. 

(22) 
(23) 



where we have used the notation /x = cos 9. Notice that, b^ 2/ means acting b on a spin-2 function 
2/ first, then acting with another b on the resulting spin-1 function b2/; the final result is a spin-0 
function, which is invariant under a rotation. The same applies to the case f' -2/- 

One can relate the spin-s spherical harmonics the usual spherical harmonics by (Goldberg et 
al. 1966) 



loi < s < I, and 



sYlm 



il + s)\\ 



1/2 



rYim 



(-l)^b-^li. 



.(^-^)!J 

for — Z < s < 0. Using these, we obtain (e.g. Hu &; White 1997b) 



(24) 
(25) 



±2Ylm 





1/2 r 


[{i+m 





2i 

di — cot OBq =b (80 — cot 

smd 



1 



sm 



Yi. 



Im- 



For completeness, the explicit expression for spin-s spherical harmonics is (Goldberg et al. 1966) 

sYue,(j)) = 



^ l 2l + l{l + my.{l-m)l 




s\n^\e/2) 



Att {l + sy.{l-s)\ 

l + s 
s — m + r 



(-1) 



l—s—r 



cot 



s— m+2r 



(0/2), 



(26) 
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where 

( P \ ^ P- 
\ Q J ~ (p-qV-Q^- 

are the binomial coefficients. 

Useful properties of the spin-weighted spherical harmonics are (e.g. Goldberg et al. 1966): 

UYlm = [{l-s){l + S + lt/'' s+lYlm. 
KYlm = -[{l + s){l-S + lp\-^Yirn, 
\>hYlm = -{l-s){l + S + l)sYim. (27) 

These relations also fix our sign/phase convention. Note that ZS97 choose a different sign convention 
for the spin weighted spherical harmonics, which results in the relation sY{^ = (—I)'* sYi-m- 



2.3. Thomson Scattering 

Now that we have set the mathematical scene we are ready to begin the description of the 
physics of CMB polarization. 



2.3.1. Scattering Matrix 

The process of scattering off a photon by a charged particle where there is no change in photon 
energy is called Rayleigh scattering. In particular, when the charged particle is an electron, the 
process is known as Thomson scattering. 

Imagine an electron at the origin for instance just before the epoch of recombination. This 
electron is accelerated by a incoming plane wave of radiation with wave vector ki and re-radiates 
an outgoing (scattered) wave along kg- We will call the plane spanned by ki and kg, the scattering 
plane (see Fig. 1). 

If the incoming radiation is polarized parallel to the scattering plane (which means that the 
electric field vector is in this plane), the differential cross section of Thomson scattering is (e.g., 
Rybicki k. Lightman 1979) 



du 



3aT I 



\ki ■ kg\', (28) 



where ctt is the Thomson cross section, and the solid angle = d{cos6)d(l) defined in the usual 
spherical coordinates. Note that the outgoing radiation is still polarized parallel to the scattering 
plane, but for right-angle scattering there is no outgoing radiation along kg since 
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x' 




Fig. 1. — Geometry of the Thomson scattering. Here ki and ks denote the wave vectors of the 
incident and scattered radiation, respectively. 



Next let us consider unpolarized radiation, which can be regarded as an independent superpo- 
sition of two linearly polarized waves with perpendicular axes. We can choose one wave as polarized 
in the scattering plane, the other perpendicular to the plane. Then the differential cross section, 
which is the sum of two polarized cross sections, is 



da 



(\k ■ + 1) ; 



(29) 



that is, the ratio of the intensities of the scattered light in directions parallel and perpendicular to 
the scattering plane is cos^ 9, where cos 6 = ki-kg. For right-angle scattering, = 7r/2, the scattered 
light is completely linearly polarized in the direction perpendicular to the scattering plane. 

Formally, we denote the incident radiation as I' = (ly, I^, C/', V), where /y and are the 
intensities parallel and perpendicular to the plane of scattering, respectively. The total intensity is 
/' = /y + and Q' = I'^ — Jy ^. The incident and scattered radiation is related by the scattering 
matrix (or "phase matrix", e.g., Chandrasekhar 1960) 



Stt 



/ cos2 9 



V 







COS 9 




\ 




cos^ j 



(30) 



(772,172,0,0) 



As an example, consider the case of right-angle scattering of unpolarized light, T' 
{I'/2){1, 1,0,0). Denoting the scattered radiation without a prime, we have 

I = (TtR • I' = aril' /2){cos^ 9,1,0,0)= aril' /2){0, 1, 0, 0); 

that is, the scattered light is polarized perpendicular to the scattering plane. We note that in 
the above expression the dimension of T is different from that of T' (due to the Thomson cross 



^The definition is somewhat arbitrary (up to a sign change). Here we follow that of Kosowsky (1996). 
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section). This is because we are considering a single scattering due to one electron. When one 
deals with the scatter by an ensemble of particles, the formal solution of the radiative transfer 
gives the scattered intensity with correct dimensionality (see e.g. Chandrasckhar 1960, §16-17). 
For simplicity, however, we shall ignore this inconsistency here and in the next section. 

Below it will be important to remember that if there is no V component in the incident 
radiation, Ray leigh/ Thomson scattering could not induce net circular polarization in the scattered 
light. 



2.3.2. Polarization at Last Scattering Surface 

Suppose the incident radiation field is unpolarized, Q' = U' = V = 0. Without loss of 
generality, we can choose the z-axis to be along the propagation direction of the scattered light, kg, 
and denote the incident radiation from {9', (f)') by the vector (recall the definition from last section) 
I'{9',4>') = I'q, ^,{1/2,1/2,0,0), where Ig, ^, is the total incident intensity. Let us also define the 
X- and y~axes to be parallel and perpendicular to the scattering plane. The scattered intensity 
vector, defined with respected to this coordinate choice, is then X = (3cjT/16vr)/g, ^,(cos^ 6', 1, 0, 0). 
Written explicitly, I{z) = (3crT/167r)/^, ,^,(1 + cos^ 6*'), Q(.z) = {Sar / IGtt) I'g, sin^ B' and U{z) = 0. 
Note that the notation reminds us that the quantities are evaluated on the chosen i-axis. 

The above results are for one single scattering. The total scattered intensities are obtained by 
integrating over all incident radiation. Since for every scattering event there is a unique scattering 
plane, thus coordinates x — y (although z is always the same), we need to specify a standard 
coordinate system, referring to which the integrated Q and U are calculated. In another words, 

3oT f j^/,, , __2, 



m = j^j dn'{l + cos'e')I'e',^', (31) 

Q{z) = ^ f dn'sm^e' cos2(f>' IL (32) 
167r J 

U{z) = [ dn'sin^ 9' sin 2(f)' I'f), a,,, (33) 

IGtt J 

where Eqn (12) is used, and —cf) is the azimuthal angle of every scattering plane in the standard 
coordinates. 

If we further expand the incident intensity by spherical harmonics 
Eqn (31) becomes 
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where we have used the exphcit form of Yoo{9', = l/VIn and ¥20(6', = V'57l67f(3 cos^ 9' - 1) 
in replacing the factor (1 + cos^ 6') in the integrand of Eqn (31). Similarly, because the integrand 
oi Q — iU is directly proportional to ¥22(0', (j)') = A/l5/327re^*'^' sin^ 6', we deduce 



(Q - iU)iz) = ^^a'i^J dJ^'sin^ 9'e^'^'¥im{ 

l.m 




Therefore, if there exist a nonzero quadrupole moment a'22 in the unpolarized incident radiation 
field, the total scattered radiation in z direction would be polarized (Kosowsky 1996). In other 
words, if we fix a coordinate system, e.g. celestial coordinates, place ourselves at the origin and 
look south, the polarization in the CMB we observe along this direction was generated by a 
moment in the incident radiation field of the electron which last scattered the light we observe. 

If we look along any other direction we see light that was scattered not along the z direction 
by some electron in the last scattering surface but some other direction n which points at us. This 
h makes an angle /3 with the z axis. We can expand the incoming radiation field in a frame whose 
z axis is pointing toward that direction, and relate the expansion coefficients with those in the 
unrotated frame (Kosowsky 1999). Denoting the quantities in this new/rotated frame with a tilde, 
we have {Q—iU)(z) oc a'22- To relate the coefficients of multiple expansion in the original/unrotated 
frame with those in the rotated frame, we note that 



I'(n) = r{hrot) = ^a'xy¥xy{nrot) 



x,y 



where n denotes the angle {9' = (3, (p') in the unrotated frame, and hrot denotes the same direction, 
but in the rotated frame; n^ot = R{P) n, where R is the rotation operator. Then 



'xy = J dflrot¥*y{hrot)I'{nrot) 

dhrot¥*y{R{p)n)l'{n) 
dnJ2l^^ymiR)Yx*min)l'in) 

m 



where Vy^ is the Wigner D-symbol (e.g. Sakurai 1985). Now, because of the transformation 
properties of the spherical harmonics under spatial rotations, only 02^ components of the incidental 
radiation field would contribute to a'22 which generates the polarization in the n direction. 

As an example let us consider an azimuthally symmetric incident field, i.e. only 7^ 0, then 
V^^{R) = dyjniP) and therefore a'22 = f^2o(/^)'^20 = {V^ / ^) (^20 ^^^"^ P (Varshalovich, Moskalev & 
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Khersonskii 1988). We obtain (Kosowsky 1999) 

{Q-^U){n) = —^-a',2 = ^^j^a,, sm (3. 

This means that the polarization due to an 020 component in the incident radiation field is maximal 
if we look along the equator, since (5 = it/2. In that case we are seeing the incident radiation field 
"edge on" . Notice also, that because the incident field is real, ago is real, and hence U {h) = 0. In 
this way the contribution from different a2m niodes of the incident radiation field to the Q and U 
we observe varies depending on our line of sight. 

There is a pictorial way to infer why the existence of polarization requires a quadrupole moment 
in the incident radiation field. Referring to Fig 2 we see that for isotropic incident radiation 
(monopolc), there is no polarization, by symmetry. For dipolc anisotropy the amplitude from the 
top is greatest, while that from bottom is weakest, and average on the sides. The polarization 
directions are horizontal for scattering from the top and bottom and vertical from the sides. The 
radiation from the top would produce horizontally polarized light, but the lack of radiation from 
the bottom cancels this effect. The net outgoing polarization is zero. We need to go to quadrupole 
moments to obtain polarization. In the configuration shown in the figure, we have a net outgoing 
horizontal polarization, perpendicular to the long lobes. 

Long before recombination, in thermal equilibrium, both polarization states of photons are 
equally populated; therefore the incident intensity field should not possess any polarization. In 
particular, V = 0. Since Thomson scattering can not generate circular polarization, as discussed 
in previous section, only Q and U polarization is expected to be present in the CMB. 

In this section we assumed that the incident radiation was unpolarized, which is true before 
recombination. In this case polarization is generated by a quadrupole moment in the incident 
radiation intensity. If we do allow polarization for the last few scatters before the photons begin to 
free-stream, there are two additional ways polarized emission can occur: polarized monopole and a 
polarized quadrupole. These two latter effects are subdominant to the first and we will not describe 
them in detail here. However, we do keep the polarization monopole and quadrupole terms when 
we discuss the solution of the cosmological photon Boltzmann equation in section 3.2. 

Gravitational wave perturbations can produce outgoing polarized radiation in additional ways 
which leads to additional terms coupling to the aoo component of the incident field as well as higher 
order terms. We include in our calculation of tensor perturbations in 3.3. 

3. Temperature and Polarization Power Spectra 

Having outlined how polarization is generated at the last scattering surface we now turn to the 
prediction of the statistical properties of the polarized CMB. Since the primordial perturbations are 
expected to be Gaussian to a high degree of accuracy and since linear theory is a highly accurate 
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Fig. 2. — By symmetry, a uniform (monopole) radiation field could not produce any polarization. 
Left: A dipole radiation field. The direction of the out-going radiation is out of the paper. The plus 
and minus signs indicate the intensity of the radiation. The light coming from both top and bottom 
will be polarized perpendicularly to the scattering plane, i.e. along the horizontal direction. The 
light from the sides will be polarized vertically. Since the sum of the incident light from vertical 
directions is equal to the sum of the incident light from horizontal directions, the net polarization 
is zero. Right: A quadrupole radiation field. The direction of the out-going radiation is also out 
of the paper. The intensity from both the top and the bottom is higher than that from the sides. 
Therefore the outgoing light has a net polarization in the horizontal direction. 
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approximation to the evolution of these perturbations until last scattering, the small anisotropies 
of temperature and polarization in the CMB are expected to follow Gaussian statistics. That 
smallness of the anisotropics is also good evidence for global isotropy and homogeneity of the 
cosmos. These two aspects, Gaussianity and isotropy lead us to believe that the statistics of the 
anisotropies are comprehensively summarized in the power spectra of the perturbations. These are 
therefore the central objects of investigation in all current missions observing CMB temperature 
and polarization. 

In this section we follow closely the work of ZS97 for the calculations of the temperature and 
polarization power spectra of the CMB.^ With the use of the spin-weighted spherical harmonics, 
ZS97 unified the treatments of these spectra. We first introduce the basic quantities necessary for 
the calculations of the perturbations that leave their signature on the CMB (§ 3.1), then detail the 
calculations of the spectra due to the scalar perturbations (§ 3.2) and tensor perturbations (§ 3.3). 



3.1. Definitions 

Suppose we are interested in the temperature and polarization anisotropy of CMB at a point 
{6, 4>) on the sky. The unit vector along the line of sight is called n. The Stokes parameters are 
measured with respect to the local coordinate system specified by the vectors (e^t, e^) on the tangent 
plane at the point of interest. Prom Eqn (12) we see that, under a rotation through an angle ip 
about n, a particular combination of Stokes Q{n) and U{fi) transforms as 

(Q ± iUy (n) = e^^'^ {Q±iU) (n) , (34) 

i.e., like spin-ib2 quantities. We therefore can expand them by the spin-2 spherical harmonics: 

{Q + iU){h) = ^ a2,im2Yim{n), 

l,m 

{Q-iU){h) = a-2,im-2Yim{n). (35) 

l,m 

On the other hand, the temperature anisotropy is expanded by the usual spherical harmonics 

'^(^) = (^T,lm^lm{n), (36) 

l,m 

since it is invariant under rotation (spin-0). 

With the help of the spin raising and lowering operators, together with Eqn (27), we obtain 
two spin-0 quantities: 

\>'^ {Q + iU){h) = ^a2,lm^\>2Yirnin) 
l,m 



''Some typographical errors in ZS97 have been corrected. 



-16- 



J2 a2,lm^{-Vil + m-^)lYlm) 
l,m 

(1 + 2)!" 



E 

l,m 



1/2 



{Q-iU){h) 




(37) 



From these the expansion coefficients can be found by using the orthogonahty of spin-0 &: 2 
spherical harmonics: 



0-T,lm 
0-2,1x11 



dn2YiUn) {Q + iU){n), 
(Z-2)!l^/2 



/ 



{l + 2)\ 

J dn.2Yi*^{n) iQ-iU)ih), 
.0 + 2)!. 



dnYil{n)\>' iQ + iU)in), 



dnYi*^{h)r {Q-iU){fi). 



(38) 



In Appendix A we discuss the parity of Q and U under spatial inversion, the operation that 
takes h into — n. We find that Q is parity even, U is parity odd. In the following we group together 
quantities of the same parity, and work with two spin-0 quantities E{n) and B{n) constructed from 
(Q + iU){n) and f [Q - iU){n) (ZS97): 

E{n) = ::l{\>\Q + iU){n) + f{Q-iU){n)) 



= E 



Lm 



{l-2)\ 



1/2 



B{n) = 



~^ (b^ {Q + iU){n)-f {Q-iU){h)) 



(39) 



2i 



l,m 



{l + 2)\ 



1/2 



L(/-2)! 



«S,im^«m("') 



(40) 



Lm 



where we have introduced coefficients which are linear combinations of a2,im and a_2,/m- 

o-EM = — («2,;m + a_2,im)/2, 
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{l + 2)\ 



a{E,B),lm- (41) 



{E,B),lm L(/-2)! 

The main advantage of working with E and B compared to Q and U is that they are scalar, 
spin-0 quantities. In this way the calculation of polarization spectra proceeds analogously to that 
of the temperature spectrum. While the same is true for the quantities defined in Eqn (37), other 

advantages of E and B arc that they arc real, and, more importantly, that they have distinct 
parities (see Appendix A). Under parity transformation the E pattern remains the same, while 
that of B changes sign. It is this property that reduces the mimbcr of power spectra (see below) 
that we need to calculate: cross-correlations between B and T or E vanish ^ (see § A). 

Finally, let us define the power spectra necessary to describe the statistics of CMB temperature 
and polarization maps: 

Cxi = 21 + 1 ^i^*x,lm ax,im)^ 

m 

" 2ITI ^^^^ 
m 

where X = T, E, B, and (...) denotes ensemble average. From Eqn (38) we see that 
{a*T,i'm' aT,im)= J Jdn'dfi {T*{h')T{n))Yi>m'{n')YUn). 

If correlation in temperature at different sky positions h' and n can be expressed as {T*{n')T{n)) = 
CTT{n' ■ n), i-e. only a function of the angle between the two position vectors, we can expand the 
correlation function in terms of the Legendre polynomial: 

Crrin' -71) = ^ '^^Cr^Pqin' ■ n). 

Q 

With the help of the addition theorem of spherical harmonics 

4-7r 



we obtain 



{a^j,^, aT,im) = '^^Tq J dO! Y*f.{h')Yii^i {fi') J d9.Yqk{n)Yim{n) 

q,k 



Cxq Sql'Skm'SqlSkm 

q,k 

Cti (43) 



Equivalently, one can define E{n) = Yimih) and B{h) = J2 Yim.{n), but these are almost identical 

to E and B defined above. 
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Similar calculations lead to other correlations between T, E, and B modes; we therefore conclude 

{C'T,l'm' (^E,lm) = ^Cl ^I'l Sm'm, 

{(^*B,l'm' aT,lm) = {"'*B,l'm' O.E,lm) = ^- (44) 



3.2. Scalar Perturbations 



Scalar perturbations (e.g. perturbations of the gravitational potential, see Efstathiou 1990, Ma 
& Bertschinger 1995, and Dodelson 2003 for detailed treatments and Bertschinger 2000 for a useful 
summary of structure formation) can be expanded into plane waves, each characterized by a wave 
vector k. We can decompose the temperature anisotropy seen by an observer at conformal time 
T into individual contributions from different k-modes of the scalar perturbation. There are only 
two relevant directions, k and the line of sight n. Without loss of generality we can always choose 
coordinates such that z \\ la. We denote this temperature anisotropy by A^(r, A;,/i), which by the 
above argument only depends on the angle between the photon direction h and k, n = h ■ k, the 
amplitude of the mode, k = |kj, and the conformal time r. The use of "A'" means the anisotropy 
is due to this single mode; and the superscript "s" denotes the scalar perturbation. 



Unlike the case of temperature anisotropy, for the polarization calculations, we need to specify 
a local coordinate system for every point on the sky. A "natural" choice is the local unit vectors 6$ 
and e<^ of the usual spherical coordinates, together with h as the third unit vector, as mentioned 
in § 3.1. This is a good choice because of two symmetries in the problem: in addition to the 
azimuthal symmetry, there is also reflection symmetries, that is, space inverse with respect to the 
plane containing k and n. Under space inverse, Q remains unchanged while U changes sign (see 
§ A). Therefore, under this choice of coordinates, the polarization can only be in either eg or 
direction, that is, purely Q with U = 0. Since V is also zero, the polarization tensor is diagonal in 
this frame (c.f. Eqn 16). 

In short, simple symmetry considerations suggest that only the Q component will be present 
(which we denote as Aq), and its amplitude depends only on k and r. We denote the polarization 
due to this single mode by Ap(r, k, ji) = Ag + iA^. Specifically, Ap 



Aq in this frame. 



The evolution of the anisotropies is described by the Boltzmann equation. In the synchronous 
gauge, we have (e.g., Ma k. Bertschinger 1995) 



Af, + zfc/xAf, = 

Ap + ifc/iAp = 

n(T,fe) = 



-Af. + A^ 



To + ^-Vb- 2^2(/x)n 



-Af. + -(1 



M^)n 



A: 



T2 



+ A^ 



P2-I 



(45) 



where h and h + Qr] are the trace and traceless scalar parts of the metric perturbation, Vb is the 
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baryon velocity, and the A^^ (j = 0, 2) are from the Legendre expansion 

I 

Using this expansion convention, n ■ Vb = —ivbH. The differential optical depth is defined through 

/TO 
K(r') dr', (46) 

where a(r) is the cosmic expansion factor, rig is the electron density and Xg is the ionization fraction. 
Note the time derivatives are all with respect to r. 

Wc caution that A|, defined here is the perturbation of photon brightness temperature, which 
is 1/4 of the perturbation in intensity; the latter is adopted by Bond & Efstathiou (1984) with the 
same notation. The equation (1) in Bond & Efstathiou (1984) is therefore 4 times Eqn (45). Also 
notice hs3 in Bond & Efstathiou (1984) is h + Arj used here. 

Note also that the three contributions to 11, the term which sources polarized photon brightness 
in Eqn (45), correspond to the contributions discussed in section 2.3.2. 

In Appendix B we show how to write down an integral solution to the Boltzmann equation 
(Seljak &; Zaldarriaga 1996). This method of solving the Boltzmann equation is called the line- 
of-sight technique, as a way of distinguishing it from traditional techniques such as the moment 
expansion. The result is 



A^(To,fe,M) = r dre-''=''S^{T,k), 
Jo 



S^{r,k) = e-«(,7 + a)+5(r)(^A^o + 2a + | + j + ^ 

J ^vu m\ mg{T) 
+9{t) « + 17 + 7712 + 



k 2A;2 / 4fe2 ' 
g{T) = ke-", 

AKro,fc,M) = j^' dTe-'-i'g{T)n{T,k\ (47) 

where the auxiliary variable x = k{To — r), and a = (h + Qf])/2k'^. In particular, the visibility 
function gij) encodes the ionization/recombination history of the Universe and is strongly peaked 
at the epoch of recombination. Also notice how simple the ^ dependences are in these solutions. 

So far all the calculations are done for one k mode; to obtain the total contribution from all 
modes we need to sum over their contributions. Consider the temperature anisotropy first. At the 
direction n the contribution from scalar perturbations is 

T%n) = J d3ke(k)Af.(ro,A;,/x), (48) 
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where ^(k) is a random variable describing the stochastic property of the initial density field; it 
satisfies 

r(ki)e(k2)) = P^(A;)<5(ki-k2), (49) 

where P*(fc) is the primordial power spectrum for scalar perturbations. Prom this we can calculate 
the power spectrum Cti (ZS97): 



Cti 



m 

^Y.({J dn'Yr^{h')T%h')y dnYr^{h)T%h)]) 



21 



21 



21 



21 



m 

X Aff (ro, k', /i') Af.(To, k, /i) (r (k') e(k)) 

I d'kP'ik) I dn'dnYim{n')Yi*^{h) 

X Af?(ro, A;,n • k) Af.(ro, A;,n • k) 

2 

^ I d'kP%k)Y I dnYilin)AUTo,k,f,) 
^ J d'kP%k)Y\£°drSHr,k) I dnYil{h)e-'^^ 
' dTS^iT,k)ji{x) 



= {A-Kf j k^dkP'ik) 
= {AT^fj edkP%k)[A^Tiik)]\ 



where 



AUk)= / drSI^{T,k)ji{x), 



(50) 



(51) 



and ji{x) is the spherical Bessel function of order I. Because of the configuration that ^ || k, the 
identity / dnY*^{h) e-''^'^ = [47r(2Z + 1)]V2 {-iyji{x)Smo has been used. 

Next let us consider the case for polarization. Because under a rotation of the plane perpen- 
dicular to n, the Stokes Q and U change into each other, we need to refer all calculations to a fixed 
coordinate. In other words, the contributions to polarization anisotropy from all scalar modes are 



{Q±iuyih) 



d^k^k) exp(^2i(/., ,) A|,(ro, k, ^) 



where (p^^ ^ is the angle needed to rotate each z \\k frame to the standard reference frame. This is 
a main source of complications associated with polarization spectra calculations (ZS97). Instead, 
in complete analogy to the temperature calculations (Eqn 48), we use (ZS97) 



d3k^(k)A|(To,A;,/x) 



(52) 
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to compute the power spectra, where the spin-0 quantity A~ is defined by Eqn (39): 

A'~{To,k,iJ.) = ^(b2A|,(To,fe,M) + f A|?(to,/c,/x)) 

-3 



Q /"TO 

= -j^ dTg{T)U{r,k)il+dlf{x'e-^^>^), 



(53) 



where we have used the fact that Ap^To, k, fi) = Ap(To,fc, /u) (recall that A^ = 0) and applied 
Eqns (22) & (23) to the integral solution Eqn (47). In obtaining the last expression, we have used the 
fact that g{db) [/(6) exp(-za6)] = g{db) [f{ida)exp{-iab)] = f{ida) \g{db) exp(-za6)] (where f{y) 
and g{y) are polynomials of y) to simplify the /x dependence. It is apparent from the definition 
Eqn (40) that 



thus 



A|(ro, k,fi) = ^ (b' A|,(ro, k, /x) - f A|?(ro, k, /x)) = 0, 
B%n) = J d3k^(k)A|(ro,fe,/x) = 0. 



Therefore the only polarization related power spectra for scalar perturbations that we need to 
calculate are Cei ai^d Cci- 

The calculation of Cei is now very similar to that of Cti- 
Cei = 21 + 1 ^^"'^M "-em) 

m 



21 + 1 
1 



2Z + 1 



1 



21 + 1 
1 



2Z + 1 



[-^ Yl f dn'dnd'k'd'kYUh')Yii{h) 

' m 

X A|(ro, k', /x') A|(ro, k, fi) (f (k') ^(k)) 

/ - 2)! 



I 5] f d'kP'ik) j dn' dnYi^{h')YC^{h)A'^{To,k,h' ■ k) A'~{To,k,h ■ k) 
j d'kP^{k)Y, \ j d^lYC^ih) j^' drg{T)Il{T,k){l+dlf{x^e-^^>^) 



= / £ dr g{r)Il{T,k){[l + dlf[x' ji{x)]) 

where wc have interchanged the order between angular and temporal integrations, and performed 
the angular integration to obtain the spherical Bessel function. To proceed, we use the differential 
equation that defines ji{x) 



j'l + 2j[/x + {l-l{l + l)/x-')ji = Q, 
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(where the primes denote the derivatives with respect to x) to replace the j'/ and terms by ji 
(ZS97): 

il + dlf[x'jlix)] = {l + dl){x'ji + d,i2xji+x^j'i)) 

= {l + dl){2ji + 2xjl + lil + l)ji) 

= (2 + /(/ + + 2xj'i + d^{{4 + l{l + - Aj'i - 2{x - l{l + l)/x)ji) 

= l{l + l){l{l + l)-2)ji/x'' 

= {l-l)l(l + l){l + 2)ji/x\ 



The power spectrum then becomes (ZS97) 

,{1 + 2)1 



Cei = i^^f'^^Jk'dkP^ik) 



dTg[T) n(T, k) — 2- 



= {Anf j k'dkP'{k)[/^%i{k)]\ 



where 



^ {^^i^y £drSUr,k)jiix), 



SUr,k) = 



3g{T)U{T,k) 



x" 



Following similar procedures we obtain the cross correlation power spectrum 

Cci = i47rf J k^dkP'ik) A|.; A|;. 



(54) 



(55) 



(56) 



(57) 



3.3. Tensor Perturbations 

The calculations of power spectra due to tensor perturbations (gravity waves) are basically 
the same, if somewhat more complicated. The reason is that, for each Fourier mode of the gravity 
waves, there are two independent "polarizations", often denoted by "+" and "x". Let ^'^ and 
denote independent random variables that characterize the statistics of the gravity waves, their 
linear combinations will prove to be useful (ZS97) 

e = {e-ie)/V2, 

e = {t+ie)/V2, (58) 

which satisfy 



(^i*(k)c^(k')) = {e*{^)eo^'))=p\k)d{k-k')/2, 

(e^*(k)e'(k')) = 0, 



(59) 
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where P^{k) is the primordial power spectrum for gravity waves. 

For each Fourier mode k, we choose z || k; in this frame, the temperature and polarization 
anisotropics due to the tensor perturbation are denoted as Ay (r, k, n) and Ap(r, k, n) = Aq +iAlj, 
where the A reminds us the contribution is from one k mode, and the superscript "t" means 
"tensor" . Note also, because there is no azimuthal symmetry present now, the anisotropies depend 
on azimuthal angle (p as well as n = k-h). Their explicit form are (Polnarev 1985; Kosowsky 1996) 

A*.(r,k,n) = [(l-/x2)e2^<^^i(k) + (l-M')e-2''^e'(k)]A^(r,fc,M), 

A*p*(r,k,n) = [{l + ^xfe^'t>e{^) + {^-^i?^-^'H\n^p{r,K^x), (60) 

where the new variables Ay(r, k, n) and Ap(r, fc, fi) are independent of ^ and satisfy the Boltzmann 
equation (Polnarev 1985) 

A^ + ife/xA^ 
A^ + ikiiK^p 

The solutions to these equations can be obtained by the line-of-sight integration (ZS97): 
A*.(To,k,n) = ((l-;x2)e2*^enk) + (l-//')e-2^<^e2(k)) j^J dr e-'-'^SUr,k), 
A*p(To,k,n) = ((l-^)2e2^'^^i(k) + (l + /x)2e-2^'^C'(k)) j^' dre-'^^^S^piT^k), 
A*p*(To,k,n) = ((l + /x)2e2^<^^i(k) + (l-/x)2e-2^^e'(k)) _^'''dre--''5MT,fe), (62) 

where 

Sir{T,k) = -he-'^ + g^, 

SUr,k) = -5*. (63) 

Next, to obtain spin-0 quantities we need to apply spin raising and lowering operators to the 
integral solutions (ZS97): 

b2A*p(To,k,n) = el(k)jJ|''°dr5MT,A:)(^-9^ + ^-^)^(l-/x2)(l-;,)2e2^■^e--'^] 

+e{^) dr S^ir, k) {-d, - ' [(1 - + ^if e-2^ V-^] 

= -<eHk)(l - /x^) e^*^ £ drSUr, k) {s{x) - iB{x)^ e"^^'' 



-/i-/i(A^-^'), 
-«(A*p + ^'), 

— A^n + - 
10 ™ 7 



3 6 3 
-ApQ + -Ap2 — — Ap4. 



(61) 
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TO / 2 ^ ^ 



+^2(k)y^ drS'p{T,k)[-d, + ,^'l [(l-^J?)(l-^ifB-^''^B-^^^] 

-C2(k)(l - /x^) e-^^"^ £° dTS\,{T, k) (£{x) - iB{x)) e-^^/^, 



'0 

where operators £{x) = —12 + (1 — 9^) — ^xdx and = 8a; + 2x'^ dx are introduced. In 

obtaining these expressions, we have worked out the results of the differentiations with respect to 
H, and used the trick /(6) exp(— iafe) = f{ida) exp(—iab) again. For example, the ^^(k) term of 
b2A^(ro,k,n): 

[-df^ + Y^)' [(l-/x')(l-/^)'e--^] = -(1 - (-12 -8z(l-/x)ar +(!-/.) V)e--'^ 

= -(1 - n^) (-12 - Six - 8xdx + x^{l- id^f) e"'^^ 
= - {8{x) - iB{x)^ (1 - //2)e-^^/^. 

Notice how these operators simplify the /x dependence of the expressions. Prom these we construct 
A~(To,k,n) and A|(To,k,n): 



At(ro,k,n) = -3^(b2A*,(ro,k,n) + fA*p*(To,k,n)) 



= (1 - /x^) (e2^'^^^(k) + 6-2^"^ ^^(k)) ^(x) dre-'^^^S^T, k), 



(64) 



A|(ro,k,n) = ^(b2A*,(To,k,n)-fA*;(To,k,n)) 



= -(1 - /x^) (e^^-^e'Ck) - e-2^'^e2(k)) H(x) £° dre'^^'^SUr, k). (65) 



The total contributions from all k modes are then 



X\n) = J d=^kA5,(ro,k,rx), (66) 

where X = T, E, B. The solution Aj,(ro, k, n) is already given in Eqn (62). These are our starting 
points in calculating the power spectra. 

First let us calculate the temperature spectrum: 

Cti = 2p~l^{(^T,lm^'^M) 
m 
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21 + 
1 

21 + 



21 + 



/•TO 

JO 

m 

X (1 - (1 - /x^) (^^"^^ dr e-^^"' S^(t, fc')) j^"" e"^^" S^(r, A;) 



/ d3kdJ7'dOy,^(n')y;;,(n)(l-(fc-n)2)(l-(fc-n')') 

m 



(i 



dr e 



-ixk-n' ot 



we will prove the two terms on the right hand side are equal. Let us denote the first term as tj, 
the second term tij. The following equations and expressions will be used: 

r.27r 

dcPe^^'^e-^'^'l' = 27r6m±2, 
the relation between Yi^ and associated Legendre polynomials (e.g. Jackson 1999) 



Ylm{lJ',<f>) 



21 + 1 {I - my.y^^ 



Att (Z + m)! 

the explicit form of associated Legendre polynomials 



p,„u) = (-ir(i - M^r/^A^PKM), 



and the relation between Pim{lj) and Pi-m{lj) 
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ti = 



We begin with the term tj: 
2tt 



21 + 

27r 



21 + 



/ k'dkP\k) 
-J k'^dkP\k) 



d(f) e 



dn 



2l + l{l-my. 



1/2 



27r 



(/ + 2)! 



k^dkP\k) 



r 



dre-^^'^5^(r, k) 



the second term t// is: 



21 + 
27r 

2rr 



k'^dkP\k) 



2r 



E 



e' 



i{—2—m)ip 



d^ 



21 + 1 {I- m)\ 
Att {l + my. 



1/2 



Pi^(/x) (1 - fi^) r dre-'^^ 5^(t, k) 
Jo 



rTX) 



dii{l-i?) Pi-2(m) j dTe-'^^'S'T{r,k) 

Pz2(/-2)!/{«+2)! 



ti = 2r 



{1-2)1 



k^ dkP\k) 



l"ro rl 

I dTStr{T,k) / dli{l- l?)Pl2{tl)^ 
J-1 



-IXfJ, 



{l + 2)\ 

Using the definition of associated Legendre polynomials, the angular integration can be written as 



d,i{l-,i^)Pi2{ii)e-'-^ = 



dji 



dfi'^ 



Pz(^) (l-MYe 



' - d^^il + din-ixe-^^>^) 
-1 J-1 



dfj, 



d/j, 



+ Pi{„){l + dinixe-'^nL 

1 

- j\f,Piii,){l + dl)\xh-^^n 
= 0+ Pi{n)i{4iijL-x + ^i^x)e-'''^'{l- 

- (1 + ^2)2 (^x^ J' diiPiiix) e-''^!'^ 



.2\2|1 



-2{l + dlf{x\-iyjl{x)l 
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where we have used 



j'^dfie-'-''Pi{li) = 2{-iyji{x). 



Therefore the temperature power spectrum is (ZS97): 
Cti = "^ti 



4t 



{l-2)\ 



dTStr{T,k)2{l + dl)\x\-ifji{x)) 



= {inf JedkP\k)[Airi{k)]^ 



(67) 



where wc have used the expression Eqn (54), the same trick used in calculation of scalar polarization 
power spectrum, and defined 



L(^ + 2)!J 



1/2 



Jo 



drS^T^r, k) 



Ji{x) 



(68) 



where S^{t, k) is given in Eqn (63). Notice the similarities between these expressions and Eqns (50) 
& (51). 

This formulation has a big payoff, because the calculations of remaining spectra are very similar 
and straightforward. The reason is that the angular dependence of A~,(ro,k, n) and A~(ro,k, n) 
are the same as that of A|,(to, k, n). This is clearly shown in Eqns (62), (64) & (65). The difference, 
namely the operators £ and B, can be applied after the angular integrations are performed (ZS97). 
This leads to 



Cei = 



^° dTSUr,k)£^'^^^ 



(47r)2 J k^dkP\k) 
{Airf j dk P\k) ° dr 5^(r, k) 



x" 



-3l{x)+3'{{x) + 



2ji{x) , 4jlix) 



+ 



x^ 



X 



, (69) 



Cbi 



(47r)2 J k^dkP\k) 
(47r)2 j k^dkP\k) 



TO 



dTS^p{T, k)B 



3 AX] 



Jo 



dTS^p{T, k) 



2j[{x) + 



(70) 



These equations can be further simplified by integrating by parts the j[ and j'^ terms to 
increase computational efficiency (ZS97). The source term Sp{T,k) = —g^ (Eqn 63) makes sure 
the boundary terms at r = vanish; the boundary terms at r = tq can be ignored for / > modes 
(i.e. excluding the monopole; see the Appendix B). The j'/ term in the integrand of Cei gives 
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-jlix)ig'^ + 25* + g^)/k^, the 4jl/x term gives -4ji{x){g'^ /kx + g'^/kx + g^/x"^); the 2j[ term 
in the integrand of Cbi gives — 2_7;(x)(5'^ + g^)lk. The final expressions are (ZS97) 

Cxi = i^TrfJ k^dkP\k)[A'xiik)?, 

Cci = {4Trf J k''dkP'{k)Airi{k)A%i{k), (71) 
where X = T, E, B; A^;(fc) is given in Eqn (68), and 

/■TO 

^\E,B)lik) = / dTS\EB){T,k)jl{x), 
J 

= a{r)[ — + ^]+'2g{r)-. (72) 



4. Summary 

We have given a detailed discussion of the physics of polarization in the CMB. Like the tem- 
perature anisotropy, polarization is anisotropic over the sky. After a review of the mathematical 
and physical background (§ 2), we have re-derived the statistical properties of the temperature and 
polarization anisotropy in terms of the power spectra in the line-of-sight formalism (§ 3, see also 
Appendices A and B), which are implemented in the code CMBFAST (ZS97). We hope that this 
uniform treatment will be of value for the beginning theorist learning about polarization and the 
solution of the linearized cosmological Boltzmann equations. We provide some suggested reading 
material in the next section. 



5. Further Reading 

In this review we follow the mathematical formalism developed in Zaldarriaga &: Seljak (1997, 
ZS97) to construct polarization power spectra. Our treatment can be complemented by the article 
by Hu k, White (1997a), who provided vivid visual aids to understanding the physics of polariza- 
tion^. A recent review (Cabella &: Kamionkowski 2004) treated the same subject under the tensor 
spherical harmonics formalism (Kamionkowski, Kosowsky, k. Stebbins 1997). 

The technique presented in ZS97 was for the full-sky in a spatially fiat Universe (which is 
strongly favored by many observations, especially those from the WMAP satellite). Zaldarriaga, 



see http://background.uchicago.edu/~whu/polar/webversion/polar.html 
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Seljak, k, Bertschinger (1998) extended the integral solutions and the calculation of polarization 
field to non-flat Universe models. Hu & White (1997b); Hu et al. (1998) unified the treatments for 
anisotropies produced by scalar, vector and tensor perturbations. If this guide is to be for beginners, 
we refer our readers to these treatments for the intermediate level. A quantum mechanically oriented 
discussion can be found in Kosowsky (1996). 

For more details on simulating the CMB polarization maps and their statistical properties, we 
refer to the discussions in ZS97 (§ V) and Kamionkowski, Kosowsky, k, Stebbins (1997). Detailed 
and general discussions of polarization experiments can be found in, e.g., Seljak (1997); Hu & 
White (1997b); Zaldarriaga (1998b); Tegmark k, de Oliveira-Costa (2001); Kovac et al. (2002); 
Kogut et al. (2003); Bond et al. (2003); Bunn et al. (2003) and references therein. For more specific 
discussions on individual polarization experiments (completed or on-going), see e.g. Keating et al. 
(2001); Hedman ct al. (2001); Kovac et al. (2002); Johnson et al. (2003); Montroy et al. (2003); 
Keating et al. (2003); de Oliveira-Costa et al. (2003); Barkats et al. (2005); Farese et al. (2004); 
Bowden et al. (2004); Cortiglioni et al. (2004); Leitch et al. (2005); Readhead et al. (2004). 

General discussions of the importance of polarization in constraining or breaking degeneracies 
in cosmological parameters, or helping to discover new physics can be found in, e.g. Zaldarriaga, 
Spergel, k Seljak (1997); Kosowsky (1999); Eisenstein, Hu k Tegmark (1999); Tegmark et al. (2004) 

and the references therein. For more specific topics, such as the reionization and polarization, see 
e.g. Hu (2000); Liu et al. (2001); Haiman k Holder (2003); for weak gravitational lensing and 
polarization, see e.g. Zaldarriaga k Seljak (1998); Stompor k Efstathiou (1999); Kesden, Cooray, 
k Kamionkowski (2002); Seljak k Hirata (2004). 

For recent general review on CMB and its polarization, see, e.g. Kosowsky (2001); Hu k 
Dodelson (2002); Zaldarriaga (2003); Challinor (2004). We also recommend the treatments in 
textbooks by Peacock (1999); Liddle k Lyth (2000); Dodelson (2003) on this topic. 

Among the many useful websites that provide up-to-date information or detailed exposures on 
CMB-related topics, we refer the reader to that of Wayne Hu^ for CMB tutorials, that of Anthony 
Banday^ for a detailed listing of CMB resources, that of Max Tegmark^ for CMB data analysis, 
those of Martin Whitc^'' and Ned Wright^^ for cosmology tutorials, and Edmund Bertschinger's 
pages-^^ for excellent lecture notes on cosmology and structure formation. An online reading list 



^http://background. uchicago.edu/~whu/indcx. html 
*http: / /www. mpa-garching.mpg.dc/~banday/CMB. html 
®http: / /www. hep. upcnn.edu/~max/ 

"'"http: / /astron. bcrkelcy.edu/~mwliitc/litnilpapcrs. html 

http://www.astro.ucla.edu/~wright/intro. html 
^^http://arcturus.mit.edu/~edbert/ 
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that also provides links to many additional useful websites can be found on Martin White's page . 

Finally, for numerical work we recommend HEALPix^^, publically available mathematical 
software for the simulation and analysis of temperature and polarization maps on the sphere (Gorski, 
Hivon & Wandelt 1999; Gorski et al. 2002). 

We thank an anonymous referee for suggestions and corrections that improve the presenta- 
tion of the manuscript, and David Larson for helpful comments. YTL thanks Jimmy Snyder for 
comments on an earlier version of the manuscript. This work has been partially supported by the 
University of Illinois at Urbana-Champaign. 



A. Parity of E and B Modes 

Here we discuss the properties of Stokes Q and U and E and B (namely E and B) under parity 
transformations. The discussion is similar to that of Zaldarriaga (1998a). Let us consider the space 
inversion r^r,6^TT — 6,^^Tr + ^. In the language of basis kets used in § 2.1.2, a polarization 
state is characterized by the expectation value of the Stokes operators. If the the basis kets are 
chosen as |ei) = eg, |e2) = e^, then they are related to their counterparts in the space-inversed 
frame by |e'^) = e'g = — |ei) and {€2) = = |e2). We therefore see that the expectation values are 
Q' = Q and U' = —U: Q has even parity, while U has odd parity. 

As for the parities of E and B, it is useful to notice that 9^' = d^p, dg' = —dg. Let n and h' 
refer to the same physical direction in the original and space-inversed frames, respectively. Prom 
{Q + iU)'{h') = {Q- iU){n) we have, using Eqn (21), 

b'(Q + iUYin') 



\>'\Q + iUyin') 



Similarly, wc have 'i'^{Q + 7T/)'(r/) =\>^iQ- iU){n). 



= -m-iV)(n), 

sm \ sm J 

= (de + ^d^) sine {-UQ-iU)in)) 

smd \ smd J 

= fiQ-iU){n). 



http://astron.berkeley.edu/~mwhite/readinglist.html 
http://www.eso.org/science/healpix/ 
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Using these results and the definition of the E and B modes, 

-2E'{h') = b\Q + iU)ih') + fiQ-iU){h'), 
~2iB{h') = \>^{Q + iU){h')-f{Q-iU){h'), 

we see that 

-2E'{n') = f{Q-iU){n) + \>'^{Q + iU){h) = -2E{n), 

-2{B'{n) = f{Q-iU){h)-\)'^{Q + iU){h) = 2iB{n). 

Therefore E and E have even parities, while B and B have odd parities. 



B. Line-of-Sight Integral Solution to Boltzmann Equation 



In this section we discuss the integral solution to the Boltzmann equation (Seljak & Zaldarriaga 
1996). We first consider the polarization due to scalar perturbations. The Boltzmann equation is 



Ap + ik/iAp 

U{T,k) 



-Ap + -(i-M')n 

At2 + Apo + Ap2, 



where, for brevity, the superscript s is omitted. If we move the term —kAp on the right-hand side 
to the other side of equality, and multiply both sides with e*^'*'^^'*, we find the whole expression 
becomes ^(Ap e*'^'^^^''). (Because of the way it is defined, de~'^/dT = k,e~'^; c.f. Eqn 46.) 
Integrating over the conformal time then gives 



d 



pro 

Jo dr 



(^Ap e*^'''^-'^^ = Ap e'''^''° - Ap e-''(^=°) 



3 r° 

4 Jo 



dr e''"'^-^ k{l - ii^)U, 



where, by definition, k(to) = 0, e"''^'^"^ = 1, and we have used the fact that k(t = 0) ^ oo (i.e. the 
optical depth to the "big bang" is infinite) , the A p e'~'*('^=o) term vanishes; after dividing through 
the factor e^''^'^°, we have 

O /-To 

Ap{To,k,fi) = t(1-/x') / dre-'^^ke-'^U, 

4 Jo 

the solution (47), where x = k{To — t). 

The solution to the temperature anisotropies is obtained analogously. We have 



Ar(ro, A;,/x) 



/•To 

Jo 



1; 1 



6 



h--{h + 6ri)P2{fi) + K A^o - m - Mf^)^ 
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We can further simplify the fi dependences in the integral by taking the advantage of the e 
term: for an arbitrary function y of r and /x, we have 

/•To /"To 1 fjp-ixn 

/ dre-"-^e-''F(r,/.) = / dr — — — e"'' y(r, /x) 
Jo Jo ikfi dT 

= |e-'^(^o)y(ro, fi) - e-^'^^^o-«(o)y(0, fi) - rdre-''"'^ (e"^ Y) ] , 
ikfi [ Jo dr } 

the first term on the right-hand side is simply Y{tq), which can be regarded as a constant contribu- 
tion to the photon temperature (the "monopole"), and can be ignored since we are only interested 
in the temperature fluctuations; the second term vanishes because of the damping before the last 
scattering. Therefore, if Y{t, jj) contains integral powers of ji, we can eliminate its /i dependences 
by successive applications of the above identity. To proceed, let us group the terms on the right- 
hand side of the integral solution by powers of /i. Recall that a = {h + 6ri)/2k'^, the first two terms 
in the square brackets reduce to rj — k^fi^ a. We therefore have 

Ar(ro,/c,/x) = ^ dre"^^^ (e-'^f] + qAto + ^9^^ - ivbg - (^k'^ae'" + ^gUj 

where we have used g = ke~'^. Integrating by parts, the term proportional to /x is replaced by 
{'^b9 + V}}g)/k, while the terms proportional to /x^ give 



(ae-'^ + 2ag + ag) + ^ (gll + 2gtl + gUj 



Grouping these in derivatives with respect to g, we recover the source term Sj,{t, k) in Eqn (47). 
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